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Abstract 

We consider a non-Hermitian Hamiltonian in order to effectively describe a two-level system 
coupled to a dissipative environment. The total Hamiltonian of the model is obtained by adding 
a general anti-Hermitian part, depending on four parameters, to the Hermitian Hamiltonian of a 
tunneling two-level system. The time evolution is formulated and derived in terms of the density 
matrix of the model, different types of decays are revealed and analyzed. In particular, the popula- 
tion difference and coherence are defined and calculated analytically. We have been able to mimic 
various physical situations with different properties, such as dephasing and vanishing population 
difference. 
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1. INTRODUCTION 



Non-Hermitian Hamiltonians with complex eigenvalues are finding numerous applications 
in modern physics Quantum scattering and transport by complex potentials 

7H9|, decaying states lOj, multiphoton ionization llMl3|. optical waveguides [14 . 



resonances 



15| are all examples of the application of non-Hermitian quantum mechanics (NHQM). It is 



also known that NHQM can be used in the theory o: 



open quantum systems, one example 



is given by the Feshbach projection formalism [10, |17|, which describes a system with a 
discrete number of energy levels interacting with a continuum of energy levels [ljj . 

In view of all these applications, the dynamics of non-Hermitian quantum systems can be 
regarded as subject of active current research (among alternative approaches to dissipation 



one could mention those based on the nonlinear evolution equations 19|, |20|). It has already 



been studied by means of the Schrodinger equation 



12 



21 



22[ or derived methods 23] 



while a different approach has been followed in [24j , where the Wigner function of an initial 
Gaussian state has been calculated considering correction terms up to quartic order in the 
Planck constant. 

It is known that a generic open quantum system can be conveniently described by a den- 
sity matrix while wave functions are to be interpreted stochastically j^l . In the present work, 
we consider the density matrix of a two-level system (TLS) and study its non-Hermitian time 
evolution in order to mimic the coupling to a dissipative environment. The total Hamilto- 
nian of the model has been obtained by adding a general anti-Hermitian part, depending 
on four parameters, to the Hermitian Hamiltonian of a tunneling two-level system. The 
dynamics of the density matrix is solved analytically in a number of relevant cases. The 
parameters are chosen in order to make the evolution free from singularity and to impose 
specific constraints on the nature of the solution. 

In section [2] we give an outline of the NHQM formalism. In section [3] we consider a 
generic two-level system, solve the evolution equations and obtain expressions for selected 
observables. In sections |4j [5] and |6] we consider special cases which correspond to different 
physical situations. Conclusions are drawn in section [7J 
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2. NON-HERMITIAN DYNAMICS 



Unless stated otherwise, the results of this section apply to a general non- Hermit ian 
system. The Hamiltonian of such system can always be partitioned into Hermitian and 
anti-Hermitian parts 

H = H + + H^ (1) 
where we denoted H± = ±H ] ± = \{H±H^). 



2.1. Evolution equations 



Upon introducing the density matrix as p(t) = \^f(t))(^f(t)\, the non-Hermitian 
Schrodinger equation, 

^|*(t)) = ~(tf + + ff_)|*(t)> > (2) 

leads to the evolution of the statistical operator in terms of a commutator and an anticom- 
mutator: 

|««) = -i [*♦,«*)] -s {*-.««)}, (3) 

where the square brackets denote the commutator and the curly brackets denote the an- 
ticommutator, respectively. As shown in j^], Eq. ([3]) can also be written in matrix form 

as 



d_ 

Of 



(4) 



where we have defined a matrix super-operator 



A 



1 + H4H+)" 1 

-1 + {H+^H- 



(5) 



and the column vector 



PA, 



Pit) 



(6) 



together with the corresponding row vector n T . . 

Since in NHQM the dynamics is not unitary, the trace of the density operator is not 
preserved in general: 

»tr(fl = ltr(A*_), (7) 



3 



hence, the operator is not a projector. Following the standard definition of averages 

for pure states in non-Hermitian quantum mechanics 1|, we introduce a normalized density 
operator 

p' = p/tr (p) . (8) 

In terms of the density operator (JSJ) the quantum average of an observable x = x(0) can be 
defined also for mixed states as 

( X )t = tr (p'(t)x(O)) = tr (p(t)x(O)) /tr (p(t)) ■ (9) 

Equation fl9]) clearly reduces to the well-known rule for calculating statistical averages in 
Hermitian quantum mechanics in all cases in which tr (p(t)) = 1. In terms of the normalized 
density operator (jSJ) the purity can also be written in a habitual form: 

P(p)=tr(p' 2 )=tr(p 2 )/(trp) 2 . (10) 



2.2. Conserved quantities 

The law of change in time of the determinant of the density operator can be found using 
the evolution equation and the matrix identity IndetM = trlnM. Hence, we obtain 

d 2 

— detp(t) = — detp(i) tr#_ . (11) 

As long as we are working in the Schrodinger representation, we can easily integrate the last 
equation: 

detp(t) = detp(0)e^ ttrf , (12) 

where the self-adjoint operator V = —iH_ will be referred as the decay rate operator through- 
out the paper. Thus, the decay rate operator with a negative trace brings a system to a 
zero-determinant state at large times (in a special case of a two-dimensional Hilbert space 
this would mean that any mixed state becomes pure at large times, regardless of an initial 
condition). If the trace vanishes then the evolution does not lead to a sink in the dynamics, 



an example belonging to the class of PT-symmetric Hamiltonians 



can be found in 



29). 



Using Eqs. ()3]) and ( fit)]) we can derive the rate of evolution of the purity: 

9~ 4 ftQ3,f) 

Wt V{p) = (13) 



where we denote 

Tl(p, f) = tr(p 2 f )trp - tr(p f )tr(p 2 ). (14) 

One can see that the purity is conserved under the general non-Hermitian evolution fl3]) if the 
condition lZ(p, T) = is satisfied. Unlike conventional quantum mechanics, this condition 
is state-dependent. 

In the case of a two-dimensional Hilbert space H2, which is of interest to the present work, 
the density matrix and the other observables are obviously represented by 2 x 2 matrices. 
Hence, one finds that TZ(p, f) can be further simplified and written in a factorized form as 

K(p, f) | rank=2 = detp (2 tr(p f) - trp tr f ) , (15) 

or, using Eq. (112]) . 

^(p,f)| rank=2 = detp(O) (2tr(pf) -trptrf) e^ trf , (16) 

so that the purity (TTOj) is conserved for any state from I-L2 whose initial density matrix has 
zero determinant. In other words, if an "H 2 state is initially pure then it stays pure during 
the non-Hermitian evolution. 



2.3. Mixed states 



Mixed state is a statistical ensemble of several pure states which are described by the 
density matrices pW = |^W)^W|. The density matrix of a mixed state can be thus written 
as the following linear combination 



P 



where the coefficients must satisfy the normalization condition 

i 

which generalizes the one used in the conventional quantum mechanics. Therefore, one can 
introduce the time- dependent functions 
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and one can also assume that pi = p'i{0) provided the corresponding density matrices have 
equal traces at t — 0. Then the normalization condition takes the habitual form 

!>*(*) = 1 ' ( 2 °) 

i 

and equation ([17]) can be rewritten as 

p'(t) = Y,P' l (t)p { Ht), (21) 

i 

where the primed operators are defined according to (JHJ). 



3. NON-HERMITIAN TWO-LEVEL SYSTEM 



We consider a system with ground and excited states denoted by \g) and \e), respectively. 



In terms of the system state projectors the Pauli operators take the standard form 



30| 



a x = \e}(g\ + \g}(e\= | ° * | , (22) 

6- y = i(\e)(g\-\g)(e\) = I ° ~' h (23) 

i 

a z = \e)(e\-\g)(g\= ( J ^ | . (24) 

The identity matrix J is given by the completeness relation in the two-state space / = 
|e) (e| + \g) (g\. In this paper we will be mainly interested in such observables as the population 
difference 

= ' (26) 

and the coherence 

<*.)t = (26) 
where Pij(t) are the zjth components of the density matrix. One can check that during the 
evolution the spin averages obey the following identity 

t~ v2 , \2 , v2 i , detp detp(O) 2 ttrf 

(^x) t + K)t + (cr*) t = 1 - 4 ^-^ = 1 - 4 ( tr ^ 2 e " < !> (27) 

which means that for pure states the averages lie on the Bloch sphere (&x)t J i~(&y)t Jr (&z)t = 1- 



We assume that in absence of any interaction with the environment (closed-system dy- 
namics), the two- level system is free to make transitions between its two energy levels. Such 
a situation is modeled by the Hermitian Hamiltonian 

H+ = -Mla x . (28) 

In order to formulate the open system dynamics of the model, we introduce a general anti- 
Hermitian Hamiltonian of the form 

H- = iMl (a I + a\b x + a 2 (Ty + a 3 a z J (29) 

and add it to the Hermitian operator defined in Eq. fl2"5]) . Equations (1251) and ( l2"9~j) can 
be substituted into Eq. ([T]) in order to obtain a total non-Hermitian Hamiltonian. It is 
convenient to choose the coefficient ao so that the density matrix is finite at all times. The 
choice of the other coefficients a^, i — 1, ...,3, is made to impose desired constraints on the 
non-Hermitian dynamics at long time. To this end, the coefficients of the ant i- Hermitian 
Hamiltonian in Eq. ( 129]) can be rewritten in a parametrized form as 

a = -7, at = 7/3, a 3 = W = yf (1 + 7 2 )(1 - P 2 ) ~ < (30) 

where the square root in W is defined up to a sign. By this we express ao in terms of 
other parameters but otherwise this reparametrization is general. The new parameters 
will be more convenient to use in what follows since a solution written in their terms has 
more physical clarity and conciseness than when using the original parameters a,. The Oj 
coefficients (i = 0, ...,3) are assumed to be real-valued. This implies that the additional 
condition 

(l+ 7 2 )(l-/3 2 )^ (31) 

is fulfilled. 

The non-Hermitian Hamiltonian studied in this work is given below 
H = H+ + iT = -Ml 

where the coefficients in Eq. ( 130]) have been used. The decay rate operator operator thus is 
given by 

t = Ml (7/3<t x - a 2 a y - Wa z - 7/) . (33) 



u x + i( 7/ - 7/3a x + a 2 a y + Wa z 



(32) 
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Considering the Hamiltonian ( 132|) and the initial condition 

P(0) = |e)(e|= I 1 °], (34) 




the corresponding equation of motion ([3]) can be solved analytically. The closed form for 
the matrix elements of the density matrix is given by 

1 e~ Tt 

= n -5* 9 [^i cos ( ut ) + B ± sin + C i cosh ( rt ) + D i sinh ( r *)] > ( 35 ) 

2 p z + 7^ 

1 e" r * 

Pi2(*) = 77-55 o [A 2 cos (wt) + S 2 sin (wt) + C 2 cosh (rt) + L> 2 sinh (rt)] , (36) 

2 p A + 7^ 

1 e" r * 

P21 (t) = - 02 , 2 [A 3 cos (wt) + B 3 sin (wt) + C 3 cosh (rt) + D 3 sinh (rt)] , (37) 

2 p z + 7^ 

1 e" r * 

P22(*) = 77755 5 [A4 cos (wt) + S 4 sin (wt) + C 4 cosh (rt) + D 4 sinh (rt)] , (38) 

2 p A + 7^ 



where 



T = 2 7 fi > 0, w = (39) 



are, respectively, the decay rate coefficient and the tunneling frequency multiplied by the /3 
parameter, and also we have introduced the following coefficients 

A x = /3 2 + 7 2 - W 2 , B x = 2/3H 7 , Ci = /3 2 + 7 2 + IV 2 , L>i = 27 W, 

A 2 = -C 2 = iW(l - a 2 + z'7/3), fi 2 = -(7 + - a 2 + z 7 /3), D 2 = iB 2 , 

A 3 = -C 3 = A* 2 , B 3 = B* 2 , D 3 = D* 2 , 

A 4 = -C 4 = -(1 - a 2 ) 2 - 7 2 /3 2 , B 4 = D 4 = 0. 

It is easy to check that the density matrix is still Hermitian but its trace is not conserved any- 
more. Such non-conservation of the trace represents the fact that the Hermitian (sub) system 
is coupled to environment, which is mimicked by the anti- Hermitian part of the Hamiltonian, 
so that the probability can be lost or gained. 

If we introduce the coefficients v4 5 = a 2 — 1 + j3 2 and C5 = 1 — a 2 +7 2 , the exact evolution 
of the trace of density matrix is given by 

tr (/W) = (*), (4°) 

where we denoted T(t) = A 5 cos (cut) + \B\ sin (ut) + C5 cosh (rt) + |Di sinh (rt). 



8 



One can check that the determinant of the density matrix given by our solution vanishes 
at all times which confirms the general formula f fl2|) since the initial density operator flMj) has 
zero determinant. If one also computes the purity (11 Op one obtains that it equals to one at 
all times, V(p(t)) = 1, which means that the non-Hermitian dynamics in this case preserves 
the purity of the initial state of the Hermitian (sub) system. It can be easily explained by 
looking at the formulae ( fl3|) . (fT6l) and (1341) . 

Using the solution for density operator and definition fl9]), one can find also the exact 
time dependence of the averages of the Pauli operators 

(a x ) t = -J— [A 6 (cos (cut) - cosh (Ft)) - 7A5 sin (ut) + (3C 5 sinh (Ft)] , (41) 
T(t) 

(e y ) t = J— [W{a 2 - 1) (cos (out) - cosh (Ft)) + (3C b sin (out) - ~/A 5 sinh (Ft)] , (42) 
T(t) 

(a z ) t = -J— [At cos (ut) + \Bi sin (ut) + C 7 cosh (Ft) + \D X sinh (Ft)] , (43) 
T(t) 

where A 6 = — r yf3W, A 7 = C 5 — W 2 and C7 = A 5 + W 2 . One can also determine the 
asymptotic value (t — > 00) of Eqs. (I41ti43p : 

lim (a x ) t = -(hny 1 lim (H+) t = 0, 

t—t+oo t— >+oo 

(a 2 -l)( 1 + W )+ 1 p 2 

hm {a y ) t = -, — . , (44) 

t^+oo" y/ a 2 -l -7(7 + W) 

,. > a 2 + lW 
hm (a z ) t = — - . 

t->+oo 7 2 + 1 

As mentioned above, the average of a z expresses the difference in the probabilities of 
occupation of the two levels while the average of a x is related to the phase of linear su- 
perpositions of the two states. Hence, the behaviour at large times of the two-level system 
with Hamiltonian defined in Eq. ( 1321) . whose time dependent properties stem from the non- 
Hermitian evolution of the density matrix, can be controlled by the parameters entering the 
definition of the decay rate operator in (1331) . In Fig.[T]we plot the analytical solutions of (a z ) t , 
(H + ) t /hQ, and (f ) t /m as functions of 2fiQt for /3 = 0.9, a 2 = 0.01 and 7//? = 0.2, 0.5, 1. 

Further, from the constraint in Eq. ( 13TT) it follows that, once 7 and a 2 are fixed, the 
oscillatory frequency is always bound from above by the critical value ou c . This latter is, 
in turn, bound from above by u + = 2Q, i.e., ou ^ uj c ^ ou + , with u c = j3 c u + , j3 2 = 
1 — (a 2 ,/^ 2 + 1)) ^ /3 2 , and j3 = ±1. This can be summarized by saying that at a given 
decay rate, the parameter a 2 measures the difference of the system's critical frequency from 
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the tunneling frequency of the Hermitian two-level system whose Hamiltonian is given in 
d28D. 



4. EVOLUTION WITH CONSERVED AVERAGE ENERGY 

There are instances in which the coupling to the environment produces dissipation while 
leaving the average energy (H+)t of the system constant. One such case is provided, for 
example, by the canonical ensemble. In order to describe the relaxation toward the constant 
average energy condition at large times, we impose 

lim (& x ) t = -(hn)' 1 lim (H+) t = 0, (45) 

t— >+oo t— >+oo 

which is equivalent, according to ()44p . to the constraint 

P = 0. (46) 

It turns out that this condition leads to the coherence (T261) in this case vanishes not just 
asymptotically but identically, 

(a x ) t = -(hQ)- 1 (H + ) t = 0, (47) 

which can be shown by directly substituting (1461) into (l4"Tj) . This model can be further cast 
into subclasses, depending on whether the parameter a 2 is chosen to be zero or not. 



4.1. Exponential decay of (cr y )t and (a z )t 

By this we assume that the observables approach their asymptotic values exponentially 
fast. If we impose that a 2 ^ then non-Hermitian Hamiltonian is defined by the sum of the 
Hermitian Hamiltonian in Eq. (I28p and of the anti-Hermitian decay operator 

f( ed ) = -Ml (a 2 a y + a 2 + 7/) • (48) 

multiplied by the imaginary unit. In such a case, the evolution equation yield the following 
solution 

p n (t) = \ [a 2 cosh (at/2) - sinh (at/2)} 2 e~ Tt , (49) 

P\2(t) = -f>2i(t) = 7^2 («2 - 1) [1 - cosh (at) + a 2 sinh (at)} e~ vt , (50) 

p 22 (t) = \ (a 2 - l) 2 sinh 2 (at/2)e~ Tt , (51) 
a 2 
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where a = 2a 2 Q and the decay rate coefficient T was defined in (|39|) . If we define 

T (ed) (t) = (al - a 2 + 1) cosh (at) - a 2 sinh (at) + a 2 - 1 , (52) 
the trace of density operator can be written as 

tr(p(t)) = le" r 'T( cd )(t). (53) 



The analytical expressions of the averages of the spin operators do not depend on 7 (yet, their 
behavior depends on the sign of 7 which determines whether the exponents are decreasing 
or increasing with time); they are given below 

(&*)t = , (54) 
( & y)t = T ( cd )( t ) (° 2 - x ) t cosn - ° 2 sinh ( at ) ~ x l ' ( 55 ) 
(& z ) t = j^ja 2 (a 2 - 1 + e~ at ) . (56) 

The asymptotic behaviour of the above-mentioned quantities significantly depends on 
values of the parameters. If, for definiteness, one assumes throughout this section that 

7 > 0, (57) 

then one can find that the parametric space of a 2 has the following physically admissible 
domains: 

(i) < a 2 < 7: 

In this case the density matrix vanishes at large times. The asymptotic values of other 
observables become: 

lim (a y ) t = -1, lim (a z ) t = 0. (58) 

t— >+oc t— >+oo 

(ii) < -a 2 < 7: 

In this case the density matrix also vanishes at large times, asymptotic values of other 
observables are given by: 

hm (cr y )t = 2 - , hm (a z ) t = (59) 
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4.2. Polynomial decay of (a y )t and (<J z )t 



By this we assume that the observables approach their asymptotic values in a way which 
is described by a polynomial or rational function of time. This can be achieved by choosing 
the decay operator as 

f(P d ) = -m (a z + 7/) . (60) 

The corresponding non-Hermitian Hamiltonian is obtained by adding ( |28|) to this decay 
operator (being multiplied by the imaginary unit). Keeping the initial density matrix in the 
form of Eq. ( 134")) . the evolution equation is solved analytically. The matrix elements have 
the following explicit time dependence 

p n (t) = (nt-l) 2 e- rt , (61) 

Pu(t) = -P21W = (^ - 1) e" r * , (62) 
p 22 (t) = Q 2 t 2 e~ rt . (63) 

If we introduce T^ pd ^(t) = 2Vtt(Vtt — 1) + 1, the trace of density operator is given by 

tr(p(t)) = e - r *T (pd) (t) . (64) 

The analytical expressions of the averages of the Pauli operators are given in this case by 
1- \ 1- \_ 2fit(l-fit) l-2fit 

and one can see that these observables evolve according to a polynomial law rather than 
exponential. Assuming that 7 > 0, the limit for t — > 00 of the density matrix and of the 
averages (cr y ) t and (& z )t are —1 and 0, respectively. This is similar to the case a 2 > above, 
with the only difference being that these observables approach their asymptotic states not 
exponentially fast but in a polynomial time. 

5. ASYMPTOTICALLY VANISHING POPULATION DIFFERENCE 

For open systems one often expects that the population difference goes to zero at large 
times: 

lim (a z ) t = 0. (66) 

t— >-+oo 
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This is equivalent, according to (|4*4|) . to the constraints 

^-7^1-/32 = 0, W = v 7 !-/? 2 , (67) 

such that the anti-Hermitian part of the Hamiltonian ( )32|) simplifies to 

f ^ = fr 7 ft [/35- a - v 7 ! - /? 2 (^ + 7"^) - i] , (68) 

and the analytical expressions for a solution from the section [3] hold provided one makes 
changes in constants according to the constraints (157)1 : 
1 e~ r * 



Pn(0 = 7T 
Pi 2 (t) 



2 /3 2 + 7 2 

1 e~ rt 

2 /3 2 + 7 2 



[A 8 cos(wt) +7 2 cosh(rt) - 2W / (/3sin(u;t) +7sinh(n))] , (69) 
[A,(cos (cut) - cosh (H)) + £ 5 (sinh (Ft) - i sin (cut))] , (70) 

(71) 



P22W = [(2 7 ^-7 2 )(cos(u;t) -cosh(n))] 

where we denoted 



7 = Vl + 7 2 , (72) 
and A 8 = 7 2 - 2W 2 , A 9 = W(t0 + i^W - 1)), 5 5 = f3(^ 2 - 7W) + 27^(7 - W)) with VP 
being defined in (I67|) . 

The asymptotic values of observables become: 



lim (a x ) t = -(HQ)- 1 lim (H + (t)) = (3, 

t— ¥+00 t— >+oo 



lim {a y ) t 

t—t+oo 



(73) 



whereas the critical frequency saturates the upper bound, uj c = u + and /3 2 = 1. We thus 
obtain 



(74) 



the profiles of observables for this model are shown in Fig. |2j 



6. DEPHASING 

It can be interesting to investigate whether non-Hermitian dynamics can cause an ini- 
tial non-diagonal density matrix to become diagonal at large times: this is known as the 
dephasing. To this end, let us choose at time zero the following density matrix 

1 v-^ / 1/2 1/2 , 

p(o) = 2 E l fc >H= ' I ■ 

k,m=g,e V 1/2 
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FIG. 2: Profiles of the averages: (a z ) (solid line), (M7) 1 (H + )t (dashed line) and (HQ) 1 (T)t 
(dotted line) versus time variable 2(3Q,t, evaluated at ft = 0.9, constraint ([67]) and, from top to 
bottom, 7//3 = 1, 1/2, 1/5. 
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If the anti-Hermitian part of the Hamiltonian ( 132]) is chosen as 



f( dph ) = htt 



o-y - 1 [cr z + I 



then the evolution equation yields 



Pn(t) 



1 



-2Tt 



Mt) = Mt)* = -e~ 2Vt [( 7 -z)e r * + z] 

p 22 w=^[(^-i) 2 +7 2 ; 



where the decay rate coefficient was defined in (139]) . 
Further, the averages can be written as 

, 2 



r 



7 2 cosh (Tt) — 1 



7 (1 



7 2 cosh (Tt) — 1 ' 

1 - cosh (Tt) - 7 2 sinh (Tt) 



7 2 cosh (Tt) — 1 

where 7 was defined in ( T72l) . We thus obtain the following asymptotics 

lim p(t) = \ (1 + 7 ~ 2 ) |<7><y|, 
lim (<7 x ) t = lim (a y ) t = 0, 



lim (a z ) t 

t— >+oo 



t—t+oo 

-1. 



(76) 

(77) 
(78) 

(79) 



(80) 
(81) 
(82) 



(83) 
(84) 
(85) 



The dephasing is accompanied by the certainty of occupying the ground state for t — > 00. 
The evolution of the observables is shown in Fig. [31 



7. CONCLUSION 

In this paper we have considered a two-level system described by a Hermitian Hamiltonian 
and added an anti-Hermitian term to it. Hence, we have considered the dynamics resulting 
from the total non- Hermitian Hamiltonian. The anti-Hermitian part of the Hamiltonian has 
been assumed to effectively describe the averaged influence of degrees of freedom associated 
with the environment. This influence is encoded in the parameters of the anti-Hermitian 
part which can be tuned in order to implement desired properties into a model. We have 
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FIG. 3: Profiles of the averages: (a z ) (solid line), (MTj^iH^t (dashed line) and (^) _1 (f) t 
(dotted line) versus time variable 2^£lt, evaluated at, from top to bottom, 7 = 4, 1, 1/5. 
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also established that the value of the trace of ant i- Hermit ian part plays an important role 
in determining the decaying behaviour of the system (or absence thereof) . 

When analyzing the evolution of the total system, we have focused on the density matrix 
of the model and expressed observables through it. Analytical solutions have been obtained 
in a number of relevant cases - such as when the evolution takes place with dephasing 
or vanishing population difference. The various types of decays produced show that non- 
Hermitian dynamics is able to mimic the influence of an environment on a two-level system. 
One of the possible future directions would be to apply this formalism to those physical 
systems which allow the two-level mode approximation. 
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